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We examine the characteristic quantities of pion-emitting sources extracted by model-independent 
imaging analysis in relativistic heavy ion collisions. The moments of the spatial separation of pion 
pair emission can provide the characteristic information about the source geometry and coherence. 
They are better for describing the non-Gaussian sources with granular and core-halo structures. An 
improved granular source model of quark-gluon plasma droplets can reproduce the main character- 
istics of the two-pion correlation functions and source functions in the experiment of ^snn = 200 
GeV Au+Au collisions. The transverse-momentum dependence of the normalized first-order mo- 
ments of the separation for the granular source is consistent with that of the usual interferometry 
results of source radii, after taking into account the Lorentz contraction in the direction of transverse 
momentum of pion pair. 

PACS numbers: 25.75.-q, 25.75.Gz 



I. INTRODUCTION 

Hanbury-Brown-Twiss (HBT) interferometry has been 
extensively used to extract the space-time and coherence 
information about the particle-emitting sources produced 
in relativistic heavy ion collisions [J 0, S 13 ■ conven- 
tional two-pion HBT analysis one needs fitting the mea- 
sured correlation functions with parameterized formulae, 
e.g., a Gaussian form, to obtain quantitatively the source 
radii and chaotic parameter. So these quantitative HBT 
results are model depended. Recently, many studies in- 
dicate that the particle-emitting sources produced in rel- 
ativistic heavy ion collisions are far from Gaussian dis- 
tributed [i i S i, 1 0, EH, E [H, Q m, [H, 113 ■ For 

the non-Gaussian sources the conventional HBT method 
of the Gaussian fit is inappropriate [H, [H, [l^ H^, HH . 
Therefore, it is important to investigate the source char- 
acteristic quantities by model- independent analysis. 

The imaging technique introduced by Brown and 
Danielewicz [a 113: HSl is a model-independent way to 
obtain the two-pion source function S{r), the proba- 
bility for emitting a pion pair with spatial separation 
r in the pair center-of-mass system (PCMS), from the 
measured two-pion correlations. This technique has 
been developed and used in analyzing one- and multi- 
dimension source geometry in relativistic heavy ion colli- 
sions 15,, ,6|, 0, a i III HH, [li, HI [2s, 2^ 21, 2a, 13, s mi . 

Recent imaging analyses for the pion-emitting sources 
produced in relativistic heavy ion collisions suggest that 
the sources may have a core-halo [gI, V\ or granular struc- 
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ture |12l | . In this paper we will examine the characteristic 
quantities extracted by the imaging analysis for the non- 
Gaussian sources with granular and core- halo structures. 
We find that the zero-order moment of r, correspond- 
ing to the intercept of the two-pion correlation func- 
tion at zero relative momentum, is usually larger than 
the chaotic parameter obtained by the HBT Gaussian 
fit. The first-order moment and standard deviation of r, 
which are normalized to the Gaussian radius for a Gaus- 
sian source, can be used to characterize the source size 
and the deviation of the source distribution from Gaus- 
sian form. 

In Ref. ^16i] , a granular source model was used to ex- 
plain the HBT puzzle, i?out/Ji!sidc '^-^ 1, in the Relativistic 
Heavy Ion CoUider (RHIC) experiments [H S, [SI . 
Here i?out and i?sidc are the HBT radii in the directions 
parallel and perpendicular to the transverse momentum 
of pion pair l36l [stI . However, in the previous granular 
source modelsflGl l38| , all the droplets are assumed with 
the same initial radius and pions are emitted at a fixed 
freeze-out temperature Tf. In this paper we will improve 
the granular source model by introducing random initial 
radii of the droplets according to a Gaussian distribu- 
tion and letting the pions emit in a wide Tf region for 
including both directly produced and decayed pions. We 
will investigate the characteristic quantities extracted by 
the three-dimension imaging analysis ^3Q . i31j | for the im- 
proved granular source. We find that the improved gran- 
ular source mode can reproduce the main characteristics 
of the experimental two-pion correlation functions and 
source functions in ^snn = 200 GeV Au-|-Au collisions 
(Tol [H, [11]. The transverse- momentum dependence of 
the normalized first-order moments agrees with that of 
the HBT radii obtained by usual Gaussian fit, after tak- 
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ing into account the Lorentz contraction in the direction 
of transverse momentum of pion pair. 

The paper is organized as follows. In section II, we dis- 
cuss briefly the limitation of the HBT Gaussian fit for the 
non-Gaussian sources with granular and core-halo struc- 
tures. The model-independent characteristic quantities 
extracted by imaging analysis are examined for the non- 
Gaussian sources and compared with the HBT results of 
the Gaussian fit. In section III we introduce an improved 
granular source model of QGP droplets to simulate the 
two-pion HBT correlations in relativistic heavy ion col- 
lisions. In section IV wc investigate the three-dimension 
source functions for the improved granular source. The 
transverse-momentum dependence of the characteristic 
quantities of the source functions are investigated in dif- 
ferent directions. Finally, a summary and conclusion is 
given in section V. 



II. GAUSSIAN FIT AND IMAGING ANALYSIS 
IN HBT INTERFEROMETRY 

In this section we discuss the quantitative results 
extracted by the usual Gaussian fit and a model- 
independent imaging analysis in HBT interferometry. 
For a direct comparison with analytical results, we con- 
sider in this section only the static sources with spherical 
symmetry. 



A. Gaussian fit 

In conventional HBT analysis, the measured two-pion 
correlation functions are fitted with a Gaussian for- 
mula to obtain quantitative results. Using the relative- 
momentum q of the pair as a variable, the Gaussian fit- 
ting formula is 



c{q) = 1 + n{q) = 1 + Ahbt eM-q^Rlht) 



(1) 



where i?hbt and Ahbt are called the HBT radius and 
chaotic parameter of the pion-emitting source. 

Let us examine the fitting procedure. The chi-square 
function is UM 
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where TZi is the ith measured correlator at qi, ct,' = (Ji/TZi, 
and CTi is the error of TZi ■ Minimizing the with respect 
to the fitting parameters i?hbt ^^'^ 
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where the quantities 

Qi = lif/cTi , b, = ln7?.i . 

Introduce 
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equations (j4]) and ([5]) can be expressed as 
1 
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It can be seen that the term for the contribution of the 
ith measured correlator has a weig ht factor of Ufqf/af. 
At a smaller qi, TZi I and di is usually larger than that 
for a larger qi . So the contributions from the data at the 
bins near zero q are very small. 

Recent imaging analyses in relativistic heavy ion colli- 
sions exhibit a two-tiered structure or a long-tail in the 
two-pion source functions 0, H, [Hlj Ell- These 

non-Gaussian source functions indicate that the pion- 
emitting sources may have a granular [T^ or core-halo 
[1, 01 structure. Let us examine the Gaussian-fit to the 
HBT correlation functions of these kinds of non-Gaussian 
sources. 

For the granular source, the pious are assumed to be 
emitted from separated droplets in the source [l^, [l^, [H, 
[39! , 40] . Assuming that the distributions of the emitting- 
points in each droplet and the droplet centers in the 
source are Gaussian forms with standard deviations a 
and i?gi- respectively, the two-pion correlation function 
for a static granular source is given by [s^, H^] 



C^gr('Z) = l + ^e- 



-q a 



-q [a +n 



(11) 



where N is the number of droplet in the source. In Eq. 
PT|) the second term is the contribution corresponding to 
the two pious from the same droplet and the third term 
is for the pious from different droplets. 

For the core-halo source l4l| , the pious are emitted 
from a central core and a halo of long-lived resonance 
decays. Assuming the emitting-points in the core and 
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halo are all the Gaussian distribution with standard de- 
viations Rc and i?h, the two-pion correlation function of 
the static core-halo source is given by jll] : 



+2/c(l - 



(1 - /c )e 



(12) 



where fc is the fraction of core emission. The second, 
third, and fourth terms in Eq. are the contributions 
for both the two pions from the core, the halo, and one 
from the core and another from the halo, respectively. 
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FIG. 1: (Color online) (a) Two-pion correlation functions for 
completely chaotic Gaussian, granular, and core-halo sources, 
(b) and (c) The fitted results of i?hbt and Ahbt of the Gaus- 
sian fit with the data sets beginning from the ith measured 
correlator. 

In Fig. [ija) we show the simulated two-pion corre- 
lations for completely chaotic Gaussian, granular, and 
core-halo sources, where the lines are the fitted curves 
of the Gaussian fit, Eq. U). It can be seen that the 
fit is very good for the Gaussian source. However, for 
the granular and core-halo sources the fits are inappro- 
priate, which miss the data of the lower g-bins. In Fig. 
[T]Jb) and (c), we exhibit the fitted results of i?hbt and 
Ahbt for the fits which use the data sets beginning from 
the ith measured correlator. One can see that the fitted 
results of i?hbt are insensitive to the first several data 
of g-bins. For the core-halo source, the i?hbt results are 
almost the same as that of the Gaussian source. The fit- 
ted HBT radii cannot reflect the large spatial extent of 
the halo for the long-lived resonance decays. Also, be- 
cause of missing the data of the first one or two g-bins, 
the results of Ahbt for the core- halo source are about 
which is the main reason that a large halo-structure was 
introduced to explain the low chaotic parameters mea- 
sured in experiments [H, |4l[. For the granular source, 
the two-tiered structure of the correlations [H, H^l and 
the failure for the first several data in the Gaussian fit 
lead to i?hbt < ^gr and Ahbt < 1- So, the usual model- 
dependent HBT fitted results are hardly to reflect the full 
geometry and coherence of these non-Gaussian sources. 



B. Imaging analysis 

The imagin g t echnique introduced by Brown and 
Danielewicz [sl |22|. |23| allows one to obtain the two-pion 
source function S{r) in PCMS, from the measured two- 
pion correlation function C{Q) [Q = \/ — Qq). After 
knowing the source function one can calculate the mo- 
ments of r model-independently. Define the nth-order 
moment of r as: 



47r drr'' S{r)r'^ 
47r dr S{r) 



1,2,... 



(13) 



where the denominator is a quantity of normalization 
as the zero-orde£ moment of r. Denote the zero-order 
moment of r by A, 



A = 47r / drS{rY ■ 
'o 



(14) 



Theoretically, the value of A is equal to TZ{Q = 0) It 
is unit for a completely chaotic source and between zero 
and unit for a partially coherent source. 

For a completely chaotic Gaussian source, the source 
function is i5„ 22x, 2:hl 
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where i?ga is called as the Gaussian radius of the source. 
The first- and second-order moments of r for the Gaus- 
sian source are: 



:i?ga , (H)ga = 6Ri, . (16) 



The deviation of the Gaussian source function is 

a ^ ^Jy^)~W = v/6 - 16/7ri?,ga • (17) 

The source functions for comple tely chaotic granular 
and core-halo sources are given by [ij, l42| 
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In Eq. (jlSp the first term is the contribution correspond- 
ing to the two pions from the same droplet and the second 



4 



term is for the two pions from different droplets, fn Eq. 
(I19p the first and second terms are the contributions for 
the two particles both from the core and from the halo, 
and the third term corresponds to that one particle from 
the core and another from the halo. 

We introduce a normalized first-order moment and a 
normalized deviation as 



R 



and 



a — a I — IG/tt . 



(20) 



(21) 



Both the quantities are normalized to the Gaussian ra- 
dius i?ga for a Gaussian source. The normalized moment 
R describes the source size. It is agrees with the HBT 
radius for a Gaussian source and we shall see that it is 
greater than the HBT radii for the granular or core-halo 
source. The normalized deviation a is equal to i? for 
a Gaussian source and larger than i? for the granular 
or core-halo source. One can examine the deviation of 
source distribution from Gaussian distribution by com- 
paring the values of R and tr. 
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FIG. 2: The two-pion source functions for the Gaussian, gran- 
ular, and core-halo sources as the same in figure [T] 

In Fig. [5] we show the source functions obtained by 
imaging analysis for the Gaussian, granular, and core- 
halo sources as in Fig. [TJ The values of A, i?, and a are 
calculated with the imaging source functions. They are 
consistent with the corresponding analytical results ob- 
tained by Eqs. ([13, (HH]) - (EH). It can be seen that the 
source function of the granular source is wider than that 



TABLE I: The characteristic quantities of the granular 
sources (J?gr = 5.5 fm). 



source parameters 


C 


i?(fm) 


5(fm) 


a = 1.2 fm, iV = 12 


9.63 


5.32 ± 0.13 


6.04 ± 0.16 


a = 1.7fm, N =12 


3.39 


5.47 ± 0.12 


6.05 ± 0.15 


a = 1.7fm, iV = 20 


1.88 


5.61 ± 0.13 


5.91 ± 0.17 



of the Gaussian source. Because of the long-lived res- 
onance halo the source function of the core-halo source 
has a very long tail. By comparing with the results of 
the Gaussian fits in Fig. [I] one can see that the results 
of A and R for the Gaussian source are consistent with 
the corresponding HBT fitted results of A^bt and i?hbtj re- 
spectively. The value of a for the Gaussian source is equal 
to that of R as expected. For the non-Gaussian sources, 
the results of A calculated with the source functions are 
close to unit and give more realistic chaotic degrees of the 
sources than that given by the Ahbt results. The results of 
R for the non-Gaussian sources are larger than the corre- 
sponding i?hbt results. The first-order moment R reflects 
more about the whole sizes of these non-Gaussian sources 
than i?hbt, because the Gaussian fit is insensitive to the 
measured data in smaller relative momentum region (cor- 
responding to larger coordinate region). The results of 
a > R for the non-Gaussian sources reflect the deviations 
of the source distributions from Gaussian form. 

The values of R and a are related to the space- 
structure of the source. For the granular source the quan- 
tity £_ = [{Rgy/a)^/{N-2)] {KN < i?^^/ of) character- 
izes the granularity of the source sensitively . In Table 
|T]we list the values of ^ as well as the results of R and 
a calculated with the imaging source functions for the 
granular sources with different source parameters. The 
values of R and a change with the droplet radius a and 
number N for the sources with the same radius Rg^. The 
difference between the values of R and a increases with 
^. The model-independent observables R and a provide 
the spatial information of the granular sources. 



III. AN IMPROVED GRANULAR SOURCE 
MODEL 

As it is seen in section II, the source functions of the 
granular and core-halo sources are wider than that of 
the Gaussian source, which agrees with recent imaging 
analyses of RHIC experiments d, 0, [13] ■ Because granu- 
lar source model can reproduce the experimental results 
of HBT radii |16j , we construct here an evolution pion- 
emitting source based on granular source model and take 
into account the effect of resonance decay by letting pions 
emit in a wide temperature region. 

On the basis of Bjorken picture [i^ the systems pro- 
duced in relativistic heavy ion collisions reach a local- 
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equilibrium at tq ~ 1 fm/c, then expand hydrodynami- 
cally. The expanding velocity of a fluid cell in a central 
rapidity region \y\ < Vm at coordinate point {r± = p, z, t) 
can be expressed as (4J] 



vi_ = v^{z = 0) Vl 



z/t. 



(22) 



We assume that the system fragments and forms a 
granular source of many QGP droplets at a time to(> 
To = 1 fm/c). The fragmentation may be due to the vi- 
olent expansion of the system with large fluctuation of 
initial matter distribution or the rapidly increased 
bulk viscosity in the QGP near the phase transition [4^ . 
Because of the surface tension of the QGP, the droplet 
has a spherical geometry in its local frame. We assume 
that the initial radii r'^ of the droplets in droplet local 
frame have a Gaussian distribution with standard devi- 
ation a, and the initial droplet centers are distributed 
within a short cylinder along the beam direction (z di- 
rection) with the probabilities [l^ 

"^^^ [l-exp(-p2/A7^i)]^?(7^^~po), (23) 



27rpo rfpo 



dyo 



|yo|), 



iotanhyo, 



(24) 



where po and zq are the initial transverse and longitudinal 
coordinates of the droplet, j/q is the initial rapidity of 
the droplet, and TZ± and ATZ± are the initial transverse 
radius and shell parameter of the granular source. In our 
calculations we take TZ± = 6.8 fm, ATZ± = 3.3 fm, and 
a — 2.3 fm. The initial temperature of the droplets is 
taken bo be 200 MeV and the quantities ym and to are 
taken to be 1.5 and 5.0 fm/c. 

The velocity of a droplet dependents on its initial 
central-coordinate (po, zq). Based on Eq. we assume 
that the velocity of the droplet is given by 



VdA. = ar 



(—) 



1 



Vdz 



zo/to 



(25) 



where ot and &t are the magnitude and exponential 
power parameters determined by the particle transverse 
momentum distribution. 

The evolution of the system after to is the superpo- 
sition of all the evolutions of the individual droplets, 
each of them is described by relativistic hydrodynamics 
with the equation of state (EOS) of the entropy density 
[46i. i47. 48]. The values of the EOS parameters we used 
are the transition temperature Tc — 170 MeV, the transi- 
tion temperature width ATc — 0.05Tc, the bag constant 
B = 250 MeV/fm^, and the ratio of the degrees of free- 
dom of the QGP to the hadronic gas dg/dn = 3.6. 

In order to include the pions emitted directly at 
hadronization and decayed from resonances later, we let 
the pions freeze-out (emit out from the source) within a 
wide temperature region with the probability 

^ (X /dh.exph(T-T^)/ATdi,.] 

+ (i-/dir)exp[-(r-r„)/Ardec], (26) 



where /dir is a fraction parameter for the direct emission, 
Tfi represents the temperature of complete hadronization, 
ATdir and ATdcc describe the widths of temperature for 
the direct and decayed pion emissions. 
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FIG. 3: The space-time distribution of pion-emitting points 
in the center-of-mass frame of droplet. 

Figure [3] shows the space-time distribution of pion- 
emitting points in the center-of-mass frame of droplet. 
Considering that most pions are emitted directly from 
the hadronization configuration and the decayed pions 
are produced within a larger temperature region, we take 
the parameters in Eq. (|26|l to be fdh- = 0.85, ATdir = 10 
MeV, and ATdoc = 90 MeV in our calculations. The 
complete hadronization temperature Th is taken to be 

0.9Tr. 
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FIG. 4: (Color online) The pion transverse momentum dis- 
tribution of the granular source and the PHENIX data for 
a/snn = 200 GeV Au+ Au collisions with minimum bias [49| . 
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FIG. 5: (Color online) (a) The two-pion correlation functions 
of the granular source and PHENIX experiment [s^l ■ (b) The 
imaging results corresponding to the correlation functions in 
(a). 



In Fig. m we exhibit the pion transverse momentum 
distribution of the granular source. By comparing with 
the experimental data [49l| , we determine the velocity pa- 
rameters of the droplet in Eq. ([25|l as ar = 0.36 and 
br — 1.50. In Figure[5](a) we show the two-pion correla- 
tion functions C{Q) of the granular source (open circle) 
and PHENIX experiment = 200 GeV Au-|-Au col- 

lisions (solid circle) [50] . In Fig. [5] (b) we show the source 
functions S{r) corresponding to the correlation functions 
in Fig. [5] (a). It can be seen that the source function 
of the granular source exhibits a long tail as the exper- 
imental result. The symbols □ in Fig. [5] (a) are for 
the restored correlation function of the granular source. 
It is consistent with the original correlation function of 
the granular source. The difference between the granular 
and experimental correlation functions at small Q may 
be arising from the final state interaction and the coher- 
ence of pion emission which we have not considered in 
our granular source model. 



IV. SOURCE HBT CHARACTERISTIC 
QUANTITIES 

In this section we examine the three-dimension imag- 
ing of the improved granular source with the technique 
proposed by Danielewicz and Pratt [13, [3l|. We will 
investigate the characteristic quantities A, R, and a in 



different directions as well as the usual HBT radii i?out, 
i?sidc, and i?iong for the granular source. 

Using a Cartesian harmonic basis {-4^ ... (il)} {I — 
0,1,2, ai = X, y, or z) [sol . Isij . the three-dimension 
correlation function TZ{Q ) = C(Q ) — 1 and source func- 
tion 5(r) in PCMS can be expressed as [l3,[33,[3l| 

7^(Q)= ^ <,....,(QX,..,o,(f^Q), (27) 



5(r) = 



E 



(28) 



where Q and r are the relative momentum and coordinate 
of the pion pair in PCMS, Q and fig are the modulus and 
solid angle of Q, and r and fir are the modulus and solid 
angle of r. From the Koonin-Pratt formalism [s^, l5l| . 
the relation of i?^ ... a,iQ) and S'^^ ... ^^(r) is given by 



,,(Q) = 4^ I drr^Ki{Q,r)Si 



,Ar), (29) 



where \M M 



ICiiQ,r) = 



0, 



for even I , 
for odd I , 



(30) 



if neglecting the final state interaction of the pion pair. 

From Eq. P7|) one can get R^oi.- - .aiiQ) from the 
measured three-dimension correlation function 7?.(Q ) 

[13, sm, 



(Q) = / ^< (^^Q)7^(Q ). (3i) 



Then, using the one-dimension imaging technique [^, 
[Hi one can get S'^^ ... „^(r) by Eq. and finally get 

the three-dimension source function S{t) with Eq. (^5)1 . 

In Fig. [5] we show the independent R^ai, -- ,ai{Q) to 
Z = 4 for the granular source. Here R^{Q) is the cor- 
relation function averaged over angles, the index "ai2". 



"ai4", and "ai2a,-2" denote "ai,ai", "a.^ 



and 



Qi, Qj", aj" , respectively. Because of source symme- 
try, the components with odd / are zero. We use x and 
y denote the directions parallel and perpendicular to the 
transverse momentum of the pion pair, and let z-axis 
along the longitudinal (beam) direction. The x, y, and z 
directions are usually called as "out" , "side" , and "long" 
directions in HBT interferometry [S^, [33| . 

The source function along the 7-axis (r — Vj, j = 
X, y, OT z) can be expressed as (lOl [30l ISlf 

S{r, ) = ^0 (r, ) + S^, (r, ) + Sf^{r,) + ■ ■ ■ , (32) 

where is the angle-averaged source function and 
S'j J (/ = 2, 4, • • • ) are the /-order modifications for S*" 
in the j-direction. 
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FIG. 6: The independent . , a, (Q) up to / = 4 for the 

granular source. 
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In Fig. [7] we show the three-dimension source func- 
tions of the granular source for kx < 350 MeV/c and 
kx > 350 MeV/c, where kr = |pit + P2t|/2 is the trans- 
verse momentum of pion pair in the longitudinally co- 
moving system (LCMS) 4]. One can see that the widths 
of the source functions in "out" direction are larger than 
those in "side" direction as observed in experiments ■ 
This reflects that the source expansion, which boosts the 
pair velocity (momentum), leads to different geometries 
in "out" and "side" directions. For higher fcy, the width 
of the source function in "out" direction is larger than 
that for lower kx- However, in "side" direction the width 
of the source function are smaller for higher kx than that 
for lower kx- The source function in longitudinal direc- 
tion has a long tail for the pion pairs with lower kx for 
the granular source. It is because that the average lon- 
gitudinal velocity of pion pairs is larger for smaller kx- 
Our model calculations indicate that the average veloc- 
ity of the pair in longitudinal direction is about 0.43 for 
kx < 350 MeV/c and 0.35 for kx > 350 MeV/c. In Fig. 
[5] we show the angle-averaged source functions S^{r) and 
the "second-order" and "fourth-order" source functions 
S^{r) and S'^(r) in "out", "side", and "long" directions 
for kx < 350 MeV/c and kx > 350 MeV/c for the granu- 
lar source. It can be seen that the one-dimension imaging 
S'" of the granular source exhibits a "two-tiered" struc- 
ture [T^. The fourth-order source functions are almost 
zero exception for that at small r. 

In order to examine the three-dimension source func- 
tions quantitatively, we introduce the moments of for 



FIG. 7: (Color online) The three-dimension source functions 
of the granular source for kr < 350 MeV/c and kr > 350 
MeV/c. rj = r^, Vy, or r^. 



the "i-component" of the source function. 



1.2, 



(33) 



and define the normalized first-order moment and devia- 
tion as 



/TT 



and 



(7,/V2-4/^, 



(34) 



(35) 



which are normalized to the Gaussian radius i?ga for a 
one-dimension Gaussian source [S{ri) ~ exp(r^/4i?ga)]. 

In Fig. [9] (a), (b), and (c) we exhibit the normalized 
first-order moments (symbols a) of the granular source 
as functions of kx in "out" , "side" , and "long" directions, 
respectively. Rx increases and Ry decreases with kx as 
expected. Because the average longitudinal velocity of 
the pairs with larger kx is smaller, Rz decreases with kx- 
Figure [D (d) gives the results of A calculated with the 
angle-averaged source function S°{r) in Eq. p4| for the 
granular source. The results of A decrease with kx- In 
Fig. [5] (e) we show the ratios '5j/Rj {j = x, y, z) for the 
granular source. The large value of Oz/ Rz at small kx in- 
dicates that the source function have a serious deviation 
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FIG. 8: (Color online) The zero-order, second-orde, and 
fourth-order source functions S'^{r), S^{r), and S'^{r) in 
"out", "side", and "long" directions for kr < 350 MeV/c 
and kr > 350 MeV/c for the granular source. 
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from that of a Gaussian source in longitudinal direction. 
However, in "side" direction the ratios are almost unit. 

The first-order moments Rx, Ry, and Rz describe the 
average separations of the source in PCMS. In LCMS the 
pair has a transverse velocity vt {— kT/[{Ei + E2)/2]). 
The spatial separation in LCMS in "out" direction is 
smaller than R^ by the Lorentz contraction factor = 
a/1 - wf, [S^]. In Fig. [S](a), the symbols v denote the 
results of ^^^Rx- One can see that j^^Rx decrease with 
kr as that of the HBT radius i?out in LCMS [l|,[3,S,l3. 
For comparison, in Fig. [9] (f) - (j) we exhibit the HBT 
radii and chaotic parameter for the granular source fitted 
with 

C((7out,fede,gio„g) = 1 + Ae-«="t«--?a-«i--?'="««'="«, 

(36) 

in LCMS and the experimental HBT results of PHENIX 
m and STAR [H for = 200 GeV Au-f Au col- 

lisions. One can see that the HBT radii as functions of 
kx for the granular source agree with the experimental 
results. Also, the separations j^^Rx, Ry, and Rz as func- 
tions of kr agree with those of the corresponding HBT 
radii i?out, ^sido, and i?iong- At lower kx the values of 
the chaotic parameter A are smaller than those of A ob- 
tained by imaging analysis, and both the values of A and 
A of the granular source are larger than the experimen- 
tal A results. The results of i?out/.Rsidc for the granular 
source are consistent with the experimental data. The 
granular source model reproduce the main characteris- 
tics of the source functions as well as HBT radii of RHIC 
experimental measurements [lol . [3^ [35| . 



FIG. 9: (Color online) (a)-(e) The characteristic quantities of 
three-dimension imaging as functions of fcr for the granular 
source, (f)-(j) The HBT Gaussian fitted results for the gran- 
ular source compared with the results of RHIC experiments 

[ll,!!!. 



V. SUMMARY AND CONCLUSION 

Ima ging analysis is a model-independent technique 
[1, [2^ |23||. With the source functions obtained by the 
three-dimension imaging technique [s^, [3l| , one can cal- 
culate numerically the first-, second-, and even higher- 
order moments of the spatial separation r in different 
directions. In principle, the detailed information about 
the source geometry, coherence, and dynamics can be ex- 
tracted by analyzing these moments. 

In this paper we examine the spatial and coherent 
information of pion-emitting source extracted by usual 
HBT Gaussian fit and imaging analysis in relativistic 
heavy ion collisions. The usual HBT results of the Gaus- 
sian fit arc model-dependent. They are inappropriate for 
describing the characteristics of the sources with non- 
Gaussian distributions, such as the sources with gran- 
ular and core-halo structures. However, the zero- and 
first-order moments as well as the deviation of two-pion 
source separation (A, R, a) obtained by imaging analysis 
can provide better descriptions for the source coherent 
and spatial characteristics. They arc model-independent 
characteristic quantities of the particle-emitting sources. 

Based on an improved granular source model we inves- 
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tigate the characteristic quantities of the pion-emitting 
source produced in relativistic heavy ion coUisions. We 
find that the granular source model of QGP droplets can 
reproduce the main characteristics of the experimental 
two-pion correlation functions and source functions in 
TiiiT = 200 GcV Au+Au colhsions 0, H HI- In 
the transverse directions of the collisions, the width of 
the source function in "out" direction is larger than that 
in "side" direction. Correspondingly, the value of the 
first-order moment is larger than that of Ry. The de- 
pendence of Rj on the transverse momentum of the pair 
kx exhibits different in "out" and "side" directions. Rx 
increases but Ry decreases with kx increase. In the lon- 
gitudinal direction, the source function for small kx has 
a long tail, which is much different from the Gaussian 
source function. Correspondingly, the value of az/Rz 
is large at small kx- We find that R^ decrease rapidly 
with kx increase. After taking into account the Lorentz 
contraction in "out" direction, all the quantities 'j^^Rx, 
Ry, and Rz decrease with kx increase. These transverse- 
momentum dependences for the granular source are con- 
sistent with those of the usual HBT Gaussian fit results 
of i?out, ^sidc, and i?iong. They are also consistent with 



the results of RHIC experiments [s^, HH] ■ 

From our model results we notice that both the values 
of A and A, obtained by the imaging analysis and HBT 
Gaussian fit, are larger than the experimental results of 
A [H, [11]. Except for the source coherence, there are 
other elements which may affect the results of A and A, 
such as the Coulomb interaction of the final-state parti- 
cles and source dynamics [li, ii,[52. Investigating 
the effects of these elements on A and A for the granular 
source will be an interesting issue. On the other hand, be- 
cause imaging analysis is performed in PCMS the tempo- 
ral information of the particle emission is hidden 5, 52]. 
Further investigation on how to extract the temporal in- 
formation of the source by analyzing the spatial moments 
in different directions and their kx dependences will be 
also of interest. 
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